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In gravity theories derived from a f(R) Lagrangian, matter is usually supposed to be 

O ' 

minimally coupled to the metric, which hence defines a "Jordan frame." However, since 
the held equations are fourth order, gravity possesses an extra degree of freedom on top 
Oh. of the standard graviton, as is manifest from its equivalent description in the conformally 
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related, Einstein, frame. We introduce explicitly this extra scalar degree of freedom in the 
action and couple it to matter, so that the original metric no longer defines a Jordan frame. 



This "detuning" puts f{R) gravity into a wider class of scalar-tensor theories. We argue 
that a "chameleon-like" detuning tracing the background matter density may provide purely 
gravitational models which account for the present acceleration of the universe and evade 



local gravity constraints. 

I. INTRODUCTION 

Observations have shown that the observable universe is well described by a Friedmann- 
Lemaitre cosmological model whose scale factor is presently accelerating [1]. The simplest way 
to account for this acceleration is to add a cosmo 



ogical constant to Einstein's equations, the 
challenge being then to explain why it is so small [2] and starts to dominate now. An even more 
economic way would be to show that this acceleration is an artefact of the averaging process, but no 
such satisfactory model has yet been proposed (see e.g. [^] for a review). As for the models giving 
a dynamical origin to this "dark energy" they fall into two broad categories: the "quintessence" 
models |^| where the acceleration is due to a scalar field coupled to gravity, and the "modified 



gravity" models [5j] where it is due to a modification of the Einstein field equations. These two 
classes of models are in fact related as most models fall into the scope of scalar-tensor theories of 
gravity jf|. 

We shall here concentrate on dark energy theories where the Lagrangian for gravity is taken 
to be a function f(R) of the Ricci scalar R [{J. As is well known [3, Q], gravity is mediated in 
such theories by the usual graviton plus a scalar, dubbed "scalaron" in [8]. However, in most 



models (see e.g. Q] for a review and recent results) matter is supposed to be minimally coupled 
to the graviton, that is, the metric, only. Following our suggestion in [lij, we propose here to 



also couple it to the scalaron, in a conformal way, so that there still exists a "Jordan frame," 
different from the original one, where matter is minimally coupled to a metric (see [ill . 



3 for 

some special examples of such non-minimal couplings). We regard this Jordan frame in which 
matter is minimally coupled to gravity as the physical frame. As we shall argue, a "chameleon-like 
detuning," tracing the background matter density, may provide a way to alleviate some problems 
faced by the minimally coupled models proposed in the literature. 

II. DETUNED f(R) GRAVITY 

Consider the following action for gravity in the presence of matter, 

S[9if,s} = ^- [ d A xV^9[f'(s)R-( S f , (s)-f(s))} + S m [^;g ij =e 2C ^9 ij ]. (1) 

Here ds 2 = gijdx l dx^ is the line element in Ai with the coordinates x l and metric gij{x l ), signature 
— h ++ and determinant g; R = g^ R%j = g^d^T^ — ■ ■ ■ is the scalar curvature, k is a coupling 
constant; the functions f(s) and C(s) of the scalaron field s{x l ) are for the moment arbitrary (we 
suppose f'(s) > 0); *(x 4 ) denotes the matter fields and S m their action. The novelty (as far as we 



are aware in this context) is that, as suggested in 10], we allow for an a priori arbitrary (conformal) 
coupling of matter to the scalaron. The field equations stemming from (JTJ) are 

f"{s)(s-R) = 2kC'(s)T, 

1 , (2) 

f'(s)Gij + -gij(sf'(s) - f(s)) + gijUf - Dijf = kT^ , 

where a prime denotes derivation with respect to s; Gij = Rij — \g%jR is the Einstein tensor; 

= -?=jlp!- is the stress-energy tensor of matter (for example, Tij = e 4C ^[(p + p)uiUj + pgij] 
for a perfect fluid, with g l3 UiUj = — 1 and p and p representing its energy density and pressure 
as measured in a locally inertial frame), and T = gijT l] . We note that T 13 is not conserved: 
D j Th = TC'{s)d l s. 

In the absence of coupling to the scalaron (C(s) = 0), we have that matter is minimally coupled 
to gij (DjT J i = 0), so that M is the "Jordan frame" representing physical spacetime, and we have 
s = R (we suppose /" / 0). Thus the equations reduce to those obtained by varying the action 
2^ J M d i x^—gf{R) + S m with respect to the metric alone. In this has been shown in fl, 

simple models such that f(R) oc 1/R n for small R may account for the present acceleration of the 
universe. 
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In presence of coupling, gravity is "detuned" in that A4 is no longer the Jordan frame. However 
JT]) can be rewritten as (we follow the notation of [13|) 



1 



S[g ij ;$] = — /_d 4 x 



M 



(3) 



where S and S are equal up to a boundary term; ds 2 = g^dxMx- 7 is the line element in the Jordan 
frame Ai with metric g~ij = Q 2C ^gij\ the potential t/(<J?) and the Brans-Dicke function cj(<3?) are 
given by (s being now a parameter) 



$(V) 



UJ(S 



3K(s)(K(s) - 2) 



U(s) = 'A')-/0O e -*7M > 



with iffs) 



dC 



(4) 



2Ws)-l) 2 " v "' d(ln/f) 

When C(s) = 0, we recover scalar-tensor theories with a vanishing Brans-Dicke function, u = 0. 

], as solar system tests impose 
We also note that in 



This property rules out models such as f(R) = R — /i 2 ( n+1 ) /R n 
u > 40, 000 



14i | when the potential can be ignored on local scales, see j^. \\\. 
regimes where K(s) — ► oo, we have lu(s) — > —3/2 and the theory resembles Palatini-/ (R) gravity 
with matter coupled to the metric only lf|. From the expression ([J]) of the Brans-Dicke function 
one may see what we are aiming at: a coupling function K(s) which vanishes on cosmological scales 
in order to account for the present-day acceleration of the universe, but which tends to 1 on local 
gravity scales in order to comply with the solar system constraints. 

Finally it will be convenient to work in the "Einstein frame" where the action (pQ) becomes 



d 4 xy/^g~* 



R* 



1 



(S» 2 -T%>) 



+ S m [* ] & j =J c Mf'(<p)- 1 gtt, (5) 
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where, again, S* differs from S by a boundary term; ds* 2 = g*jdx % dx J is the line element in the 
Einstein frame M* with metric g*j\ C(tp) and /'(</?) stand for C(s((p)) and f'(s(<p))\ the scalar field 
(p and the potential V(ip) = V(s(ip)) are given under a parametric form as 

sf'(s) - f(s) 



<p(s) = V31n v / 7 7 (i), V(s) 



4/' 2 ( S ) 



(6) 



In this Einstein frame the field equations ([2]) read 

dV 1 - K(<p) 
□ if - — = 
dip 



kT* 



G* 



The stress-energy tensor, T* 1 



2^/3 

2di<pd j( p + <7* [(a» 2 + 2V(tp)] = kT*. . 



J 3\ 
2 SS„ 



(7) 

(8) 



(for example, T£ = e 4C ^ /' M~ 2 [(/5 + p)u* iU * + pg. 



for a perfect fluid, with g*^u*u* = -1), is not conserved: L>*T*^ = -[(1 - K(<p))/VS\T*dnp. The 
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function K(<p) = K(s(<p)) is defined by (@D and ©. When K = Eqs. dSHSD are the Einstein 

frame version of the standard, minimally coupled f(R) gravity. In the "maximally detuned" case 

I I 

K = 1 (advocated in e.g. [171 )) the Jordan and Einstein frames coincide, and the theory reduces to 
General Relativity minimally coupled to a scalar field with a potential given by ([6]). For example, 
if f(R) oc 1/R n , V{ip) decays exponentially and the theory is akin to the scaling field models of 
dark energy first studied in [3]. Finally, for C(s) oc ln/'(s), that is, K = constant (see 12]), the 
theory resembles the non-minimally coupled quintessence models [191 ]. including the "chameleon" 
models [20 ], albeit with a potential which is typically exponential for large <p and hence different 



from the fiducial quintessential model 



2jJ: V(<p) = M 4 e M "^ n . 



III. COSMOLOGICAL EVOLUTION 

When the metric is taken to be that of a spatially flat Robertson-Walker spacetime, ds* 2 = 
— dt 2 + a* 2 (t)dx 2 , and matter to be pure dust, the Einstein frame field equations (j7HE|) reduce to 

Q[ r*2 -2 r rf s * ■• , o rr* ■ , dV 1 ~ K {f) * ■,, * P*0 e^"^ . . 

3H -ip - 2V(<p) =Kp , tp + 3H <p + — = — 2 ^ up with p = , (9) 

where H* = a* /a* and a dot representing derivative with respect to t. The observable Jordan 
frame scale factor is given by a(t) = a*e c( -^~^/"^ where t is the cosmic time in the Jordan frame 
(dt = e c ^-^dt). 

Suppose that, as in Q], f(s) w — l/(^ 2 s n ) (n > 0) for small s, where s = £ 2 s with t being of the 
order of the Hubble radius today. Then <p is large and V(ip) oc e~ Xlp with A = -^(n + 2)/(n + 1). 
When matter has become negligible, the Einstein frame scale factor and the scalar field do not 
depend on C((p), and they are given by a* oc t q with q = 3(n + l) 2 /(n + 2) 2 and <p ~ pint with 
p = \/3(n + l)/(n+2) [4, 5, 18]. On the other hand the Jordan frame scale factor a(t) does depend 
on the coupling function C(f). If C(ip) — > for large cp, then we have S(t) oc f( 2ra + 1 )( n + 1 )/( n + 2 ) j us t 
as in the minimally coupled case [5|, and the required accelerated expansion is achieved for n > 2. 
If now K = i^DE is taken to be a non-zero constant for large tp (that is, C(ip) ~ (^de/ V / 3) ( / 9 )) 
then the Jordan scale factor is given by 

3 V^ + ly V 2 ^ + 1 + ^DE(n + 2)y V 7 

Thus there appears another possibility, -ftToE ~ — l/( n + 1) 5 to have late time accelerated expansion 
besides JTde ~ 0, n large. 

Of course, for the cosmological scenario to be viable, this late time accelerated expansion must 
be preceded by a matter dominated era with scale factor 5 oc t 2 ^ . Suppose first that at this matter 
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dominated stage, the function f(s) is such that we can ignore V(<p) in (J9j) . This corresponds to 



the so-called "99-MDE" phase [22|. We also assume that the function K = Km is approximately 
a constant (that is, C(ip) ~ (Km/V^)<p)- Setting a* oc t q and ip ~ pint, one finds from (|9|) that 
q = 6/((l - ET M ) 2 + 9), p = 2V3(1 - # M )/((1 " ^m) 2 + 9), and 

2(1 -if M ) 2 
3(3 - (1 - K M f 



a oc t 3 ( 1+ »M) with w M = oto M ™ > 2 > • (11) 



Thus, in the "tuned" case when Km = 0, we obtain the behaviour a oc t 1 / 2 during the matter 
era, which is unacceptable 22J. The same disqualifying result holds for the detuned case if the 
matter era takes place when the potential has already reached its asymptotic regime, that is when 
s is small so that f(s) ~ — l/(£ 2 s n ) and hence <p is already large, since the condition for late 
acceleration is i^T ~ (or K ~ — l/(n + l)) in that regime. Therefore, in scenarios when the matter 



231 ]. detuning is of no help. 



era takes place when 93 has already reached its asymptotic regime 

However, if the matter era takes place when ip is small (this should be realized by an appropriate 
form of f(s) at large s (= s£ 2 )) then we may easily devise a function K(ip) which evolves from 
Km ~ 1 fo r large (p to Km ~ for small (p. In that case we have a oc i 2 / 3 during the matter phase. 
More specifically, let us take 

"W"^-^^' (12) 

so that K = 1 — (3<p + • • • at (p <C 1. At leading order, we have p = 0, hence ip is constant, (p = cpM 
(<C 1), and the Jordan and Einstein frames coincide (a ~ a* and t ~ t). Then p to first order is 
given by 

p= w (13) 

It may be noted that this result is perfectly consistent with Eq. To the accuracy of our interest, 
we have q = 2/3, hence 

3F =«p v , + 3ffv = ^_ = _. (14) 

Setting </9 = p/t, the second equation gives p which is consistent with ()13[) . 

We note that the only thing we have shown so far is the existence of two eras, the "right" matter 
dominated stage and the accelerated stage at late times. What remains to be done is to study in 
detail the dynamical system ([9]) to see which classes of functions f(s) and C(s) properly connect 
the matter era to the late accelerated era. This is left for future studies. 
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IV. LOCAL GRAVITY CONSTRAINTS 

The standard way to proceed is first to choose a background solution of the equations of mo- 
tion, ([7]) and (Jlj), then linearise and solve for the perturbations and finally check that the linear 



approximation was valid [g, [15(] . It is now well established that if we take the accelerating cosmo- 

logical solution of the model f(R) = R — /x 2 ( n+1 ) j 'R n [5( as the background, then the Eddington 

parameter is 7 = 1/2 [15| , in gross violation with the solar system observations of the Shapiro time 

delay where it was found that 7 = 1 to better than 10 -5 that is u) > 40,000. If we follow 

this standard procedure, the same result would be obtained also in the detuned case, since the 

cosmological solution is the same as in (al. 

□ 

Now, as argued in [20], the relevant background when studying the solar system is not the 
cosmological solution but the solution of (UJ) corresponding to the galactic environment. Consider 
a static, spatially homogeneous solution <p = cp g of the Klein-Gordon equation in a background of 
a uniform density p = p g . Then we have 

dV l--g(yg) * n ... * e 4C , 

Suppose (|15p has a solution for (p g small and K(tp g ) close to unity on galactic scale. For example, 
this can be realized by C(ip) given by (fl~2j) . In this case (fl~5l) gives 

Thus unless there is an extra length scale in the potential V other than the present Hubble scale, 
the right hand side is very small, typically of order p c jp g = 0(1O~ 5 ) where p c is the cosmological 
density. As for the Einstein equation ([8]), its solution on solar scale is almost flat. Thus we may 
assume the flat background. 

We now linearise the Klein-Gordon equation ([7]) around the background solution, that is we set 



(17) 



ip = (p g + ipi and g*, = f'e 2Cg (%' + hij). We obtain [6|, ll 



[A — m )ipi « j= —p Q with G eS = — and m = -^r-r-^ 

V3 87r /' e z °e dip 2 



where p~Q is the density of the Sun and A the standard 3-Laplacian. Typically (see the example 
below) m 2 = 0(np g ) or smaller. Since A = 0(L~ 2 ) 3> m 2 , where L is the typical solar system 
scale, the scalaron mass m is negligible. 

We emphasise that this is where our "chameleon detuning" differs from the standard one. 
Indeed, in the standard f(R) case the Brans-Dicke function is zero (or too small as in [12]]) and 
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one has to invoke, on top of the galactic environment, a complicated nonlinear effect inside the 
Sun. Here, to the contrary, the Brans-Dicke function oj is large and the mass m is small so that 
we remain within the regime of the linear approximation. 

For negligible to 2 , the solution of (fTTj) is ip\ « G e ff (1 — K^)Mq/ '(v3r). Since ipi <C we check 
that we could use the linear approximation. 

As for the Einstein equation (jHJ), the linearised equation SG*j ps gives the (linearised) 

Schwarzschild metric (see [g] for a detailed calculation), 



f 

ds 2 - /g 



g2C g 



(18) 



yielding the Jordan frame metric: 



with G = G e ff (1 + (1 - K g y/3) and 



d? = ^; = -(i-^)tf + (i + »^)^, ( 19) 

2 , 



2(1 - 1C)2_ 

^2 



^-jw^- (20) 

Hence for standard, minimally coupled f(R) models with K = 0, we recover 7 = 1/2 [la ]. We also 
see that, in order to comply with the Cassini mission result we must have |1 — K g \ < 0.01. 

Using "detuning" we have therefore found a way to evade local gravity constraints, using one 
aspect of the "chameleon" mechanism 



20 



23], that is the fact that the solar system galactic 



environment is much denser than the cosmological background, while remaining within the well 
understood linear approximation. 

In fact, this "chameleon- like" detuning may even help to reconcile the most emblematic model 
of f(R) dark energy Q], 

m = Us-±) so that %)~e-^ (e ^_ 1) ^ I) (21) 

with local gravity constraints. For the sake of the example, let us consider the following coupling 
function: 

c{s) = ^W thatis K{s) = J(s) with f' = e2v/v ~ 3 - W 

As can easily be seen, we have then that p* ps p and ip g ps \^3/(2£ 2 Kp~ g ) = 0(p c /p g ) for large n. 
Hence m 2 ps d 2 V7d(/? 2 | g ~ —(2n\^£ 2 ip g )~ 1 = 0(np g /n) (m 2 < but is small, yielding a Dolgov- 
Kawasaki instability [24j] setting in on a galactic time scale). Finally 1 — K g « 2<p g /V3 so that 
7 - 1 ps -8<p 2 J9 = O(10- 10 ) well below the Cassini bound. 
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V. CONCLUSION 



We have introduced a chameleon-like coupling of matter to the scalaron degree of freedom 
of f(R) gravity and argued that such non- minimal coupling might retain the late cosmological 
acceleration they predict while rendering them compatible with local gravity constraints. Of course 
many aspects of these "detuned" f(R) gravity theories remain to be explored, starting with a 



detailed analysis of the dynamical system ([9]) in the line of 
using f(R) functions such as those proposed in 251 ] . 

While we were writing up this paper, we became aware of 



231 ] and its cosmological perturbations, 



present proposal, the coupling function K (called Q in 121 ]) is a given constant 



121 ] where, in contrast with the 
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